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An innovative test for detecting quantum chaos based on the analysis of the spectral fluctuations 
regarded as a time series has been recently proposed. According to this test, the fluctuations of a 
fully chaotic system should exhibit 1// noise, whereas for an integrable system this noise should obey 
the l// 2 power law. In this letter, we show that there is a family of well-known integrable systems, 
namely spin chains of Haldane-Shastry type, whose spectral fluctuations decay instead as l// 4 . We 
present a simple theoretical justification of this fact, and propose an alternative characterization 
of quantum chaos versus integrability formulated directly in terms of the power spectrum of the 
spacings of the unfolded spectrum. 
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In the absence of a formal definition of quantum chaos, 
several different tools have been used in the literature 
to distinguish the chaotic versus integrable character 
of a quantum system. The most widely used among 
them is the study of the density of normalized spacings 
between consecutive levels of the "unfolded" spectrum, 
which gives a quantitative estimate of the local level fluc- 
tuations. More precisely, for quantum systems whose 
corresponding classical limit is chaotic, the celebrated 
Bohigas-Giannoni-Schmit conjecture |l| posits that this 
density should coincide with that of a suitable Gaussian 
ensemble in random matrix theory (RMT). On the other 
hand, a long-standing conjecture of Berry and Tabor 
states that the spacings distribution of a "generic" inte- 
grable quantum system should follow Poisson's law. Both 
of these conjectures have been verified for a large num- 
ber of systems, both numerically and using semiclassical 
analytic methods Q. 

More recently, an alternative characterization of quan- 
tum chaos was proposed in an interesting paper by Re- 
laho et al. 0]. The basic idea behind this characteriza- 
tion is considering the sequence of energy levels as a time 
series, whose corresponding power spectrum is then ana- 
lyzed using standard procedures. It was observed in the 
latter paper that for all three classical random matrix 
ensembles the (averaged) power spectrum of the fluctua- 
tions of the spacings exhibits 1// noise, while the noise 
of a spectrum with a Poissonian spacings distribution be- 
haves as l// 2 (the so-called Brown noise) . These numer- 
ical observations were theoretically explained using the 
machinery of RMT in Ref. [|, where it was also conjec- 
tured that the 1/f (resp. I// 2 ) noise of the power spec- 
trum is in fact a universal property of all fully chaotic 
(resp. integrable) quantum systems. Further numerical 
simulations have lent additional support to this conjec- 
ture. More precisely, the 1/f law has also been detected 
in the two-body random ensemble [6j], whereas a more 
general power law 1//" with 1 ^ a ^ 2 is observed 
for quantum billiards (or potentials) and random matrix 
ensembles interpolating between fully chaotic and inte- 
grable regimes, with a respectively attaining the values 
1 and 2 in these limits Q- 



An essential ingredient in the theoretical justifica- 
tion of the 1/f 2 law for the power spectrum of inte- 
grable systems is the Poissonian character of the spac- 
ings distribution I5|. However, in a series of recent pa- 
pers [U, IS US [Tllfl^l it has been shown that there a 
is wide class of integrable quantum models whose spac- 
ings distribution is not Poissonian, namely spin chains of 
Haldane-Shastry (HS) type These models, which 

are the prime examples of integrable spin chains with 
long-range interactions, appear in connection with sev- 
eral phenomena of physical interest such as strongly cor- 
related systems [TJ, generalized exclusion statistics [l5j . 
and the AdS-CFT correspondence [Hj]. It is therefore 
natural to ascertain whether the 1/f 2 law conjectured in 
Ref. @ applies to spin chains of HS type. These chains 
are particularly good candidates for a time series analy- 
sis, since their spectrum can be exactly computed for a 
very large number of sites, and there is ample numerical 
and theoretical evidence that their level density (whose 
knowledge is essential for unfolding the spectrum) be- 
comes Gaussian when the number of sites is sufficiently 
large d IS E E3, E E3, IS ■ In this letter we shall focus 
on the simplest type of HS chains, whose raw spectrum 
is equally spaced (as in the Polychronakos chain) or al- 
most equally spaced (as in the original Haldane-Shastry 
chain). Our main result is that the power spectrum of 
these chains — unlike all the integrable models previously 
studied in the literature — clearly obeys the 1/f 4 power 
law, characteristic of black noise. In order to highlight 
the main ideas and simplify the theoretical derivation we 
shall primarily concentrate on the case of an equispaced 
raw spectrum, giving only a qualitative justification of 
our result in the general case. 

Let Ei < ■ ■ ■ < E n+ i be a spectrum of a quantum 
system, and denote by \i and a its mean and standard 
deviation. In accordance with the previous remark, we 
shall assume that the following conditions are satisfied: 

(i) The energies are equispaced, i.e., Ej+i — Ex + jd. 

(ii) The continuous part of the level density is a Gaus- 
sian with parameters \x and <r, with a 3> \fnd. 

The requirement that a ^S> \fnd, which shall be needed 
in what follows, is satisfied by all spin chains of HS type 



2 



studied so far when the number of sites is sufficiently 
large [T^, El, EH- By the second condition, the unfolded 
spectrum is given by = F(Ei), where 
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The i-th normalized spacing of the unfolded spectrum is 
then defined by Sj = n{ti + i — £i)/(e„ + i — ei), so that the 
mean spacing is one. Following [J], we shall characterize 
the fluctuations of the spacings from their mean by the 



statistic 8, 



i), i = i. 



Fourier transform is given by 
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, n, whose discrete 



k = l,...,n. (1) 
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In this letter we shall be interested in the behavior of the 
power spectrum of the statistic 8i , defined by 

P(k) = \8 k \ 2 , k = l,...,n. 

Note that, since the Si's are real, 8 n -k is the complex 
conjugate of 8k, so that P(n—k) = P(k). For this reason, 
from now on we shall take k in the range 1,2,..., [n/2]. 

To begin with, we shall prove a remarkable identity 
expressing P(k) directly in terms of the Fourier transform 
of the spacings Sj. Indeed, taking into account that 
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from Eq. |T]) we obtain 
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The first term in the previous formula can be expressed 
as Y]j=i f{j)sj, with 
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Since X)J=i s i = n i the previous equality and Eq. @ 
yield the exact identity 



P(k) 
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It should be emphasized that Eq. ([3]) is in fact valid for 
any finite spectrum. The latter equation shall be the 
starting point to derive an analytic approximation to the 
power spectrum of an energy spectrum satisfying condi- 
tions (i) and (ii) above. 

In practice, to smooth out spurious fluctuations one 
usually divides the whole spectrum into several subspec- 
tra with an equal number of levels, and determines the 



average (P(k)) of the individual power spectra. For this 
reason, we shall consider a subspectrum with energies 
Ej = Ei + (j — l)d, where j = lo, . . . , l\ = Iq + m 
and to 3> 1. Calling v — m/(e; 1 — e; ), the spacings 
Si = v(ei 0+i - e/ 0+i _i), i = 1, . . . , to, are approximately 
given by 
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The discrete Fourier transform of these spacings is then 
vd 



Sk 
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When k -C to, the previous sum can be approximated 
with great accuracy by an integral, so that 
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Setting Xj — {Ei j — n)/(y/2a), y — ^/2irak/{md), Zj = 
Xj + iy, with j = 0, 1, and using [lj| Eq. 7.4.32] and 
Eq. ((3]), we finally obtain 
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(4) 



If one is dealing with the whole spectrum, then to = n, 
erf x\ ~ 1, erf a; — — 1 and Eq. (0| simplifies to 
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We shall next use Eq. Q to determine the asymptotic 
behavior of P(k) in the range fco <C fc -C to, where 
fco = (md/a) max(l, |xo|, |a;i|) <C to by the condition 
c S> TnA Since fe <C to, we can write (J4j) as 



P(fc) 
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where (p(k) = e 2y | erf zo — erf zi| 2 . The condition fc ^> 
ko implies that 

|zj|>y>l, y»|ay|, j = 0,l, (7) 

so that the well-known asymptotic formula [ijj 7.1.23] 
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holds for both z = zq and z — z\. Taking into account 
that 2y(xi — xq) = 2nk G Z we easily obtain 
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If xl 



which is only possible when x\ 



-x > 0, 



i.e., when the subspectrum under consideration is ap- 
proximately symmetric with respect to //, Eq. (|HJ) imme- 
diately yields 
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Substituting into Eq. © we finally obtain the asymptotic 
power law 
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In particular, if we are dealing with the full spectrum we 
have the slightly simpler relation 
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In the generic case, i.e., when x\ ^ x 2 , the leading term 
of the expansion of Eq. ([H]) in powers of 1/y becomes 
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Substituting again into Eq 

p(k) 
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6l) we obtain 
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When studying the whole spectrum the previous formula 
simplifies to 
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Equation (|11[) is the key analytic result for establishing 
the l// 4 behavior of the averaged power spectrum (P(k)) 
for spin chains of HS type. Indeed, when one averages 
the individual P(k) of a large number of subspectra, at 
most one of these subspectra can be symmetric about \x. 
Hence, if the total number of subspectra is sufficiently 
large we must have 
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where m is the number of spacings in each subspectrum. 

In order to test our theoretical results, we shall first 
consider the ferromagnetic rational chain of type B and 
spin 1/2 introduced in (20l |. whose Hamiltonian is given 
by 



H' 



£ 



SiSjSij 
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(14) 

Here (3 > 0, e = ±1, Sy is the operator which permutes 
the i-th and j-th spins, Si is the operator reversing the 



i-th spin, and £j = \/2ti, where t, is the i-th zero of the 
generalized Lag uerre polynomial L^~' 



_ . _ ,, 4 ...... . _/jy . It has recently 

been shown [10| that the partition function of this model 
is simply given by 
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so that condition (i) above is satisfied with 

Ei=0, d=l, n=-N(N+l), (16) 

It was rigorously proved in Ref. [HI that when the num- 
ber of sites N tends to infinity the level density tends to 
a Gaussian with parameters 

M =^(iV + l), CT 2 = ^(AT + i)(iV + l). (17) 

The spectrum is obviously symmetric about its mean, so 
that if we consider the whole spectrum P{k) should ap- 
proximately be given by for all k <C N 2 , and follow 
the power law (jTUJ) when N <C k <C A 2 . We have ver- 
ified that this prediction is in excellent agreement with 
the numerical data for a wide range of values of A, up 
to A — 100. For instance, in the latter case it is ap- 
parent from Fig. [I] that the approximation (jTDJ) is ex- 
tremely accurate even for k close to [n/2]. As to the 
I// 6 power law, the relative root mean square (RMS) 
error in the fit of the log 10 of the RHS of Eq. (TIT))) to 
log 10 P (fc) in the range 50 ^ fc ^ 1000 ~ n/5 is only 
6.9 x 10 -4 . The behavior of P(k) for small fc seen in 
Fig. [U perhaps deserves a brief explanation. To this end, 
note that when fc <C {nd/cr) rnin(|xo|, \xi\), which in this 
case is tantamount to fc <C fco = O(N), instead of ([7]) 
we have y <C min(|xo|, |xi|). Since, by Eqs. (fT6j) and 
(flT). |x | = -x and |x x | = Xi are (^(A 1 / 2 ), when A is 
large enough both — zq and zi are close to the positive 
real axis and their modules are at least 0(N 1 ^ 2 ), so that 
erf(— zq), erf z\ ~ 1. From Eq. ([5]) we then obtain 



P(k) 



-2y 



4sin 2 (^) 



fc < fc • 



(18) 



This is indeed an excellent approximation to P(k) for 
small fc, as can be seen from Fig. [TJ 

Even though the spectrum of the chain (TTJ} is symmet- 
ric about n, the l// 4 power law predicted for spin chains 
of HS type clearly emerges when one considers the av- 
eraged power spectrum (P(k)) of even a relatively small 
number of subspectra. To be more precise, we have di- 
vided the energy range [/i— 5<r, fi+5a] into ten subspectra 
of equal length m = [a] and computed numerically the 
resulting (P(k)) for a wide range of values of A, up to 
A = 500. As can be seen in Fig. [21 Eq. (TIB"]) is in excel- 
lent agreement with the numerical data except for fc close 
to [m/2] in the logarithmic scale. The accuracy of this 
approximation steadily improves as A increases; for in- 
stance, the relative RMS error of the approximation (TTj 
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FIG. 1: log-log plot of P(k) for the whole spectrum of the 
spin chain (|14[) with 100 sites (blue), compared to its analytic 
approximation ([5]) (red), the power law (|10[) (green line), and 
the approximation (|18[) (gray curve). Note that the blue and 
red dots are practically indistinguishable for log 10 k < 3. 



in a log-log plot for 10 sj k < [to/6] and N = 100, 250, 
and 500 is respectively given by 0.065, 0.046, and 0.033. 
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FIG. 2: log-log plot of (P(k)) for the spin chain lfI3|l 
(blue marks), compared to its analytic approximation (|13|) 
(red lines). Note that the initial (curved) part of the plot 
of P(k) observed in Fig. [T] is absent in this case, since 
(md/a) min(|xo|, |a;i|) < 1 for all subspectra involved. 

Consider next the spin 1/2 (antiferromagnetic) 
Haldane-Shastry chain, with Hamiltonian 



1 + 5, 



2 , * — sin" 



J ~- (19) 
A V ' 



Although the partition function of this chain has been 
evaluated in closed form Q, in practice its complicated 
structure precludes the calculation of the spectrum for 
N > 40. On the other hand, the mean and variance of 
the energy can be exactly computed for all A by taking 
suitable traces in Eq. (flSty , with the result [1] 



N 



(N 2 - 1) 



N 
480 



(A 2 - 1)(A 2 + 11) . (20) 



It was also numerically shown in the latter reference that 
for sufficiently large A the level density is well approx- 



imated by a normal distribution with the above param- 
eters. As is well known, the spectrum of the Haldane- 
Shastry chain (fT9|) is not equispaced. However, for large 
A condition (i) above is approximately satisfied, since the 
vast majority of the differences di = Ei+± — Ei are equal 
to d = 1 (for even A) or d — 2 (for odd A), and the dif- 
ferences di ^ d actually correspond to energies Ei in the 
tail of the Gaussian level density; cf. Ref. [TlJ|. As to con- 
dition (ii), the requirement a ^> \fnd is also fulfilled for 
large A, since a = 0(A 5 / 2 ), while it was known in [TlT | 
that n ^ E n+ i — E\ = 0(N 3 ). For these reasons, it is 
to be expected that the power spectrum of the Haldanc- 
Shastry chain obey the 1// power law when the number 
of sites is sufficiently large. We have numerically verified 
that this law clearly holds when N ranges from 26 to 
36. For instance, in Fig. [3] we present the results for 
TV = 32, 34, 36, where (P(k)) was computed by averaging 
over ten subspectra of equal maximal length to in the in- 
terval \fx — 5fi, [i + 5er]. A least squares fit of log 10 (P(fc)) 
to P — alog 10 fc in the range 1 ^ k < [m/6] yields an 
optimum a of 3.884, 3.908, and 3.933 for N = 32, 34, 
and 36, respectively. The corresponding values of the 
squared correlation coefficient r 2 are 0.9975, 0.9989, and 
0.9996, which strongly suggests that the l// 4 power law 
also holds in this case when A is sufficiently large. 
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FIG. 3: Plot of log 10 (P(fc)) for the spin chain dHJ) (blue 
marks) fitted by a straight line in the range ^ log 10 k ^ 
log 10 [m/6] (red). The plots for N = 32 and N = 36 have 
been displaced to avoid overlapping. 

In conclusion, the results of this letter show that 
there is a whole family of integrable systems, namely 
spin chains of Haldane-Shastry type, the fluctuations of 
whose spectrum clearly exhibit l// 4 black noise rather 
than the l// 2 noise conjectured in Ref. Note, how- 
ever, that our findings do not invalidate the theoretical 
justification for the l// 2 law proposed in the latter ref- 
erence. Indeed, an essential assumption of this justifi- 
cation is that the spacings distribution be Poissonian, 
which is certainly not the case for spin chains of HS 
type d, d, |T(| HH [3 ■ Given the fact that the spectrum 
of a chaotic system features l/f a noise with 1 ^ a < 2, 
and that spin chains of HS type possess a higher de- 
gree of integrability than generic integrable systems due 
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to their underlying Yangian symmetry, it is tempting to 
conclude that the exponent a in the 1// Q spectral noise 
provides a quantitative measure of the degree of integra- 
bility of a system. It would be of considerable interest in 
this respect to find integrable quantum systems featur- 
ing 1//" noise with 2 < a ^ 4. Finally, a noteworthy 
by-product of our analysis is the universal identity ([3]), 
which shows that P(k) oc \sk\ 2 /k 2 for k <C n. This fact 
makes it possible to translate any statement about the 
behavior of the power spectrum of the statistic Si into a 
simpler statement on the power spectrum of the spacings 
Sj. From this alternative point of view, generic integrable 
systems (with Poissonian spacings) are characterized by 



white noise in the fluctuations of their spacings s,, while 
in all chaotic systems the corresponding noise actually 
grows with the frequency. 
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